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Abstract. We consider the effect of noise on the dynamics generated by volume-preserving maps 
on a d-dimensional torus. The quantity we use to measure the irreversibility of the dynamics is the 
dissipation time. We focus on the asymptotic behaviour of this time in the limit of small noise. We 
derive universal lower and upper bounds for the dissipation time in terms of various properties of 
the map and its associated propagators: spectral properties, local expansivity, and global mixing 
properties. We show that the dissipation is slow for a general class of non-weakly-mixing maps; on 
the opposite, it is fast for a large class of exponentially mixing systems which include uniformly 
expanding maps and Anosov diffeomorphisms. 



1. Introduction 

The origin of irreversibility in dynamical systems is often modeled by small stochastic perturba- 
tions of the otherwise reversible dynamics. These perturbations may be attributed to uncontrolled 
interactions with the 'environment', or with internal variables neglected in the equations. In ex- 
perimental or numerical investigations, stochasticity or 'noise' is introduced respectively by finite 
precision of the preparation and measurement apparatus, and by rounding-off errors due to finite 
computer precision. 

One may take these stochastic perturbations explicitly into account by adding a term of Langevin 
type in the evolution equations, or equivalently introducing some diffusion term in the Fokker- 
Planck equation. In the present article we choose to deal with discrete-time dynamics on some 
compact phase space, namely the d-dimensional torus. All the maps we will study preserve the 
Lebesgue measure, which is therefore the "natural" invariant measure associated with the dynam- 
ics. The noise is implemented through a convolution operator, the kernel of which has a (small) 
width e > 0. For a given map F, two types of stochastic perturbations will be considered: the 
noise operator may be applied at each step of evolution, resulting in a "noisy evolution" ; on the 
opposite, we may choose to introduce some stochasticity only at the initial (preparation) and final 
(measurement) steps, resulting in a "coarse-grained evolution". 

In general, the influence of the noise on the long-time evolution of the system depends on the 
dynamical properties of the map. Typically, the effect is stronger if the map is 'chaotic' than if 
it is 'regular' |22IS1|2S]- We will put this statement into quantitative estimates. To this end, 
we will characterize the effect of noise on a map F through a single observable, the dissipation 
time. By dissipation, we mean the damping (through the noisy or coarse-grained evolution) of the 
density fluctuations, while the total probability remains constant. Mathematically, this damping is 
expressed by the decay of the L 2 -norm of the density (or equivalently, the L 2 -norm of the density 
fluctuations) with respect the the Lebesgue measure. The decrease of this norm due to dissipation 
may be interpreted as an "entropy increase" of the system, up to the state of maximal entropy, which 
is the uniform density (here the word "entropy" denotes the Boltzmann entropy of the probability 
density, unrelated with the topological or Kolmogorov-Sinai entropies associated with the map). 
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We define the dissipation time as the time needed for the evolution to bring fluctuations under a 
fixed threshold (that is, reduce the L 2 norm by a fixed factor) [Tl]. We are especially interested in 
the behaviour of this time in the limit of small noise. We will show that this behaviour drastically 
depends on the dynamics generated by the noiseless map F: in short, the dissipation will be "fast" 
for a chaotic dynamics, as opposed to "slow" for a regular dynamics, the difference being ever 
more striking as the level of noise decreases. Our main results exhibit this opposition by means of 
quantitative estimates on the asymptotics of the dissipation time. 

Since the dissipation time is defined in terms of an I? norm, it is naturally related to the 
spectral properties of the propagator associated with the map, acting on the space of square- 
integrable functions. In Section |3J we analyze the links between, on one side, the spectrum and 
pseudospectrum of the noisy or noiseless propagators, and on the other side, their dissipation time. 
The relevance of the pseudospectrum for time evolution problems has been recently put in evidence 
in the context of non-unitary continuous-time dynamical systems |12j . These spectral relationships 
will be mainly used to analyze the case of "regular", precisely non- weakly-mixing maps. On the 
opposite, they are of little help for more chaotic ones. 

In the following Sections, we connect the dissipation time to more 'dynamical' properties of the 
map F (under some smoothness assumptions on F). We first obtain lower bounds for the dissipation 
time from the local expansion properties of the map: a weak, or absence of local expansion will 
imply slow dissipation (Section HJ. Next, we obtain upper bounds for the dissipation time using 
information on mixing properties of the map (that is, the time decay of correlations between 
observables). Strong (e.g. exponential) mixing is a typical characteristics of chaotic behaviour and 
can be easily measured in numerical simulations [3]. The main conclusion is that exponential or 
stronger mixing implies fast dissipation, for both the noisy and the coarse-grained evolution. 

In Section El we describe several families of volume-preserving maps on the torus, for which 
the results obtained above yield useful information concerning the dissipation. The main two 
families are, on the one hand, expanding (noninvertible) maps, on the other hand, Anosov (or at 
least partially hyperbolic) diffeomorphisms of the torus. All these maps are at least exponentially 
mixing, so the dissipation is "fast" . In case of Anosov maps, we collect our results in Theorem^ As 
a matter of illustration, we also analyze in detail examples of linear mixing maps, for which exact 
asymptotics of the dissipation times can be obtained, and therefore give an idea of the sharpness 
of the bounds obtained before (the results concerning linear automorphisms had been obtained in 

EI)- 

From these developments one notices that dissipation time provides a robust characteristics of the 
chaoticity of a given dynamical system. As opposed to e.g. the decay rate of dynamical correlation, 
dissipation time has the same asymptotics (up to a constant factor) among Anosov diffeomorphisms 
while the decorrelation may be exponential (generic Anosov case) or super-exponential depending 
on particular map (cf. results regarding hyperbolic toral automorphisms |3j and Corollary EJ) . 

2. Setup and notation 

2.1. Evolution operators. Let (T d , ¥>(T d ), m) denote the d-dimensional torus, equipped with its 
o— field of Borel sets and the Lebesgue measure m. Let F : T d — > T d be a map on the torus 
preserving the Lebesgue measure: for any set B E I$(T d ) we have m(F^ 1 (B)) = m(B). In general, 
F is not supposed to be invertible. In the following we call such a map 'volume preserving' with 
implicit reference to the Lebesgue measure. 

The map F generates a discrete time dynamics on T d , which in terms of pathwise description can 
be represented by the forward trajectory {F n (xo), n E N} of any initial point (particle) xq E T d . 
However, instead of looking at the evolution of a single particle, one can consider the statistical 
description of the dynamics, that is the evolution of a density (more generally a measure) describing 
the initial statistical configuration of the system. 
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Let M{T d ) denote the set of all Borel measures on T d . For any fj, G Ai(T d ) and / G C°(T d ) we 
write 



= / f(x)d f i(x). 

The map i* 1 induces a map -F* on A4(T d ) given by 

(F»(/) = M (/ o F), for all / G C°(T d ). 

This map can also be defined as follows: 

(F*n)(B) = fiiF-^B)), for all B G B(T d ). 

In particular if \x = 5 X0 then -F*(/u) = 5f(;eo) an d one recovers the pathwise description. 

If fj, is absolutely continuous w.r.t. m, then F*(fi) preserves this property (since the measure- 
preserving map F is nonsingular w.r.t. m, see |2H| p.42]). The corresponding densities g = -j^ G 
L 1 (T d ) are transformed by the Frobenius-Perron or transfer operator Pp jH]: 



^dm J dm 
If the map F is invertible, Pp is given explicitly by: 

(P F g)(x) = (g o F'^W^p^ix) = g o F -1 (a;). 

If the map i 7 is differentiable, and the preimage set of x is finite for all x, the Perron- Frobenius 
operator is given by 

(P F g)(x)= £ 



where Jf{v) is the Jacobian of F at y. 

On the other hand one can consider the dual of the Frobenius-Perron operator, called the Koop- 
man operator, which governs the evolution of observables / G ^^(T^) instead of that of densities 
g G L . The Koopman operator Up is defined as 

(1) U F f = f° F. 

Due to the nonseparability of the Banach space L°°(T d ), it is often more convenient to consider its 
closure in some weaker LP norm, which yields larger (but separable) spaces of observables L p (T d ). 
In this paper we will be mainly concerned with the space L 2 (T d ) and its codimension-1 subspace 
of zero-mean functions Lq(T") = {/ G L 2 (T 2 ) : m(f) = 0}. This subspace is obviously invariant 
under Up and Pp, due to the assumption F*m = m. Throughout the paper, || • || will always 
refer to the L -norm (and corresponding operator norm) on Ll(T d ) (any other norm will carry an 
explicit subscript). 

For any measure-preserving map F, the operator Up is isometric on L 2 (T d ) and L^^T d ). When 
F is invertible, Uf is unitary on these spaces, and satisfies Up = Pp 1 = Pp-i. 



2.1.1. Additional notation. Although the operators introduced in previous sections were mostly 
3( 



defined on the space Lq(T^), it will be useful to consider other function spaces, which we define 



now in some detail. For any m G N, we denote by C m (T d ) the space of m-times continuously 
differentiable functions, with the norm 

||/||0m= ]T Halloo 
\ct\±<m 
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(we use the norm |oc|i = a± + . . . + for the multiindex a G N d ). For any s = m + n with 
m = [s] G N, 7] G (0, 1), let C*(T d ) denote the space of C m functions for which the m-derivatives 
are r/-H61der continuous; this space is equipped with the norm 

II /II llfll V- \D a f(x)-D a f(y)\ 
||/||c« = ll/llc- + 2^ sup —— 

\ a \ 1= m x *y lx Vl 
The Fourier transforms of functions g G L l (M. d ) and / £ L 1 (T d ) are defined as follows: 

(2) V£GlR d , 5(0 = / sC^e" 2 ™^, 

(3) VfcGZ rf , />) = / /( a; )e- 2 — fc ^ = ( efc ,/). 

Above we used the Fourier modes on the torus e^ix) := e 2nlx k . For any s > 0, we denote by 
H s (T d ) and .fP(IR rf ) the Sobolev spaces of s-times weakly differentiable L 2 -functions equipped with 
the norms || • \\h s defined respectively by 

\\9\\ 2 H . m = I (i + \a\ 2 y\m\ 2 dii, 
ii/Hh.(i-) = E( 1 + i fe i 2 ) s i/»i 2 - 

kei d 

Finally, for any of these spaces, adding the subscript will mean that we consider the ([/^-invariant) 
subspace of functions with zero average, e.g. C^T^) = {/ G C J (T d ), m(f ) = 0}. 

2.2. Noise operator. To construct the noise operator we first define the noise generating density 
i.e. an arbitrary probability density function g G L 1 (IR c( ) with even parity w.r.t. the origin: 
g(x) = g(—x). The noise width (or noise level) will be given by a single nonnegative parameter, 
which we call e. To each e > corresponds the noise kernel on M. d : 

g £ (x) = - d g (*) , 

with the convention that go = 5q. The noise kernel on the torus is obtained by periodizing g e , 
yielding the periodic kernel 

(4) g e (x) = ^ 9e(x + n). 

We remark that the Fourier transform of g e is related to that of g by the identities g e (k) = g € (k) = 
g(ek). 

The noise operator G e is defined on any function / G L^(T d ) as the convolution: 

GJ = g e *f. 

As a convolution operator defined by an L 1 density, G e is compact on -Lq(T^) (if g is square- 
integrable, G e is Hilbert-Schmidt). The Fourier modes {e k , k G Z d \ {0}} form an orthonormal 
basis of eigenvectors of G e , yielding the following spectral decomposition: 

(5) V/GL 2 (T d ), GJ= £ 9(tk)(e k ,f)e k . 

o^fcez d 

This formula shows that the eigenvalue associated with e k is g(ek). Since g is a symmetric function, 
this eigenvalue is real, so that G e is a self-adjoint operator. Its spectral radius r sp (G e ) is therefore 
given by 

(6) r sp (G e ) = \\G e \\ = sup \g(ek)\. 

4 



Since the density g is positive, g attains its maximum g(0) = 1 at the origin and nowhere else. 
Besides, because g G L 1 (M rf ), g is a continuous function vanishing at infinity. As a result, for small 
enough e > 0, the supremum on the RHS of © is reached at some point ek close to the origin, and 
this maximum is strictly smaller than 1. This shows that the operator G e is strictly contracting on 
Ll{T d ): 

(7) Ve>0, ||a|| =r sp (G e ) < 1. 

In the next section we study this noise operator more precisely, starting from appropriate assump- 
tions on the noise generating density. 

2.3. Noise kernel estimates. In this subsection we present some estimates regarding the noise 
operator, which will be used throughout the paper to estimate the dissipation time. 

We will be interested in the behaviour of the system in the limit of small noise level, that is the 
limit e — > 0. It will hence be useful to introduce the following asymptotic notation. Given two 
variables a t , b e depending on e > 0, we write 

(8) a € < b € if limsup-^ < oo, 

e->o b e 

(9) a € b € if lim = 1, 

e^o b e 

(10) a e ~ b e if a e < b t and b e < a e . 

In order to obtain interesting estimates on the noise operator G e , it will be necessary to impose 
some additional conditions on its generating density g, regarding e.g. its rate of decay at infinity, 
or the behaviour of its Fourier transform near the origin. 

The weakest condition considered in this paper is the existence of some positive moment of g, 
by which we mean that for some a G (0, 2], 

(11) M a = \x\ a g{x)dx < oo 

(we take the length |x| = (xf + . . . + x^) 1 / 2 on This condition implies the following properties 
of the Fourier transform g (proved in Appendix I A. 1J) : 

Lemma 1. For any a G (0,2] there exists a universal constant C a such that, if a normalized 
density g satisfies then the following inequalities hold: 

(12) V£ G M. d , < 1 - g(0 < C a M a \£\ a . 
If HM) holds with a = 2, we have the more precise behaviour: 

1 — ~ l£| 2 ^ n limit £, — > 0. 
In the case a < 2, we will sometimes assume a stronger property than (j!2j) . namely that 

(13) 1 - g(£) ~ \£\ a in the limit ^ -» 0. 

Notice that this behaviour implies a uniform bound 1 — g(£) < C|^| 7 for any 7 < a and C 
independent of 7. 

Typical examples of noise kernels satisfying ()13|) include the Gaussian kernel and more general 
symmetric a— stable kernels |271 p. 152] defined for a G (0,2]: 

(14) g e , a (x) := £ e-^ a/2 e k (x), 



where Q denotes an arbitrary positive definite quadratic form. For the values of a indicated, the 
function g e ,a{%) is positive on R rf . 
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In view of Eq.©, the properties or (|13j) directly constrain the rate at which G t contracts 
on LQ(T d ). For instance, (|13|) implies that in the limit e — > 0, 

(15) l-||G £ ||~e Q . 

The following proposition describes the effect of the noise on various types of observables, in the 
limit of small noise level. The proofs are given in Appendix IA.2I 

Proposition 1. i) For any noise generating density g £ ^(W 1 ) and any observable f £ LQ(T d ), 
one has 

(16) ^/-/ll^o. 

To obtain information on the speed of convergence, we need to impose constraints on both the noise 
kernel and the observable. 

ii) If for some a £ (0,2] the kernel g satisfies or then for any 7 > there exists a 

constant C > such that for any observable f E ^(T^), 

(17) \\G e f-f\\<Ce^ a \\f\\ H -r^, 

where 7 A a := mm{ / y,a}. If f E C l (T d ), the above upper bound can be replaced by 

(18) \\GJ - /|| < Ce 1Aa ||V/|| < Ce 1Aa \\Vf\U 

Using the noise operator, we are now in position to define the noisy (resp. the coarse-grained) 
dynamics generated by a measure-preserving map F. 

2.4. Noisy evolution operator and dissipation time. The noisy evolution through the map 
F is constructed by successively applying the Koopman operator Uf and the noise operator G e , 
therefore by taking powers of the noisy propagator 

T £ = G e U F . 

In general, the operator T e is not normal, but satisfies r sp (T e ) < \\T e \\ = \\G e \\. 

We will also consider a coarse-grained dynamics defined by the application of the noise kernel only 
at the beginning and the end of the evolution. Hence we define the following family of operators: 

= G e U£G e , n E N. 

In view of the contracting properties of G e , the inequalities ||T e n || < ||G e || n , HT^H < ||G e || 2 
imply that both noisy and coarse-grained operators are strictly contracting on L^{T d ). Our aim is 
to characterize the speed of contraction of these two evolutions, that is the behaviour of the norms 
||T™||, ||T e I) in the joint limits e — > and n — > 00. This characterization will be connected with 
dynamical properties of the map F. 

There are many ways to measure this speed of contraction. For instance, for fixed e > 0, the long- 
time decay of ||T e n || may be super-exponential (in which case T t is quasinilpotent) or exponential, 
governed e.g. by the largest eigenvalue of T e . However, such exponential behaviour may appear 
only after a transient time. In this paper we will characterize the noisy dynamics through a single, 
robust characteristics, namely the dissipation time. In its general form the dissipation time t* is 
defined in terms of the norm || ■ || Pj o ° n the space L^T^), and an arbitrary threshold rj E (0, 1): 

(19) Ti' v (e) := mm{n £ N : ||T e n || Pi o <v}, 1 < P < 00. 

We will be concerned with the behaviour of the dissipation time when the level of noise becomes 
small (as we prove in Proposition [3 this time diverges in this limit). In ^1] it was shown that this 
asymptotic behaviour is independent of the choice of < i] < 1 and 1 < p < 00. Therefore, in the 
present paper we will consider the computationally convenient choice p = 2 and rj = e _1 . We will 
henceforth drop the superscripts, and the dependence of on e will always be implicit: 

(20) r* := T^ e ~\e) = min{n G N : ||T e n || < e" 1 }. 
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A similar dissipation time will be defined for the coarse grained evolution: 

(21) f* := min{n G N : ||f e (n) || < e" 1 }. 

The dissipation time does not depend on whether the dynamics is applied to densities (i.e. by 
the Frobenius-Perron operator) or to observables (by the Koopman operator). Indeed, the norm of 
an operator equals the norm of its adjoint |3U1 p. 195], so that 

\\fW\\ = \\G £ U F G £ \\ = \\(G £ U F G £ )*\\ = \\G e P$G e \\, 

and similarly for the noisy operator T t . In particular, for invertible maps the dissipation time does 
not depend on the direction of time. 

We will distinguish two qualitatively different asymptotic behaviours of dissipation time in the 
limit e — > 0. We say that the operator T e (or the map F associated with it) respectively has 

I) simple or power-law dissipation time if there exists (5 > such that 

r* ~ 1/e", 

II) fast or logarithmic dissipation time if 

r* ~ ln(l/e). 

We will also talk about slow dissipation time whenever there exists some (3 > s.t. 

r* > 1/e?. 

In case of logarithmic dissipation time, the dissipation rate constant R*, when it exists, is defined 
as 

(22) H„ = lim ■ 



ln(l/e) 

A similar terminology will be applied to the coarse-grained dissipation time f*. 

3. Spectral properties and dissipation time of non-weakly-mixing maps 

In this section we investigate the connection between the dissipation time of the noisy propagator 
T e and its pseudospectrum together with some spectral properties of Up and G t . All the operators 
considered in this section are defined on L^(T d ). In the framework of continuous-time dynamics, 
connections have been obtained between, on one side, the pseudospectrum of the (non-selfadjoint) 
generator A, and on the other side, the norm of the evolution operator e tA |12| . We will obtain 
results of the same flavor, yet the proofs seem here easier than in the case of continuous time. 

3.1. Definitions and general bounds. Let us start with the definition of the pseudospectrum 
of a bounded operator |29j . 



Definition 1. Let T be a bounded linear operator on a Hilbert space Tt (we note T G C{TL)). For 
any 5 > 0, the 5 -pseudospectrum of T (denoted by a$(T)) can be defined in the following three 
equivalent ways: 

||( A _ rr i|| 

3v G H, \\v\\ = 1, \\(T-X)v\\ < 5}, 
3B G £{TC), \\B\\ < 5, \£a{T + B)}. 



(I) a s (T) = {A G C 
(//) a s (T) = {A G C 
(///) a s (T) = {\£C 



We will apply these definitions to the operator T t . For brevity, the resolvent of this operator will 
be denoted by -R e (A) = (A — T e ) _1 . We call S r the circle {A G C : |A| = r} in the complex plane, 
and define the following pseudospectrum distance function: 



d e {r) := mf{5 > : a 5 (T e ) n S r ^ 0}. 
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From the definition (I) of the pseudospectrum, one easily shows that this distance is also given by 

(23) d-\r) = sup ||i? e (A)||. 

|A|=r 

We first establish general (abstract) bounds for the dissipation time in terms of the spectral 
properties of G t and T e . In a second step, we relate these properties to dynamical properties of the 
underlying map F. 

Theorem 1. For any isometric operator U on Lq(T^) and noise operator G € , the dissipation time 
of the noisy evolution operator T e = GJJ satisfies the following estimates: 

1 - e~ l 1 

(24) - — — < r* < + 1, 

1 j 4(1) " " |ln(||G e ||)| + ' 

(25) n< inf 

r sp (T e )<r<l | ln(r) | \d e (r) 

We notice that the first upper bound does not depend on U at all, but only on the noise. Using 
the estimate (|15JI . we obtain the following obvious corollary: 

Corollary 1. // the noise generating density satisfies the estimate (|13j) for some a € (0, 2], then 
for any measure-preserving map F the noisy dissipation time is bounded from above by: 

r* < e~ a . 

Proof of the Theorem. 

1. Lower bound 

We use the following series expansion of the resolvent |3(Jl p. 2 11] valid for any |A| > r sp (T e ): 

oo 

(26) R € (\) = ^A-^Tf. 

n=0 

Considering that r sp (T e ) < \\G e \\ < 1, we may take |A| = 1, and cut this sum into two parts: 

T*-l 

R e (\) = ^ A -n-1 T" + \~ T *T[*R e (\). 

n=0 

Taking norms and applying the triangle inequality, we get 

T*-l 

\\Re(X)\\ < || A-"- 1 Tf|| + |A|--||T-||||^(A)|| 

n=0 

< n+e~ 1 ||i? e (A)|| 
=> ||i? e (A)||(l-e- 1 ) < T,. 
Taking the supremum over A € S± yields the lower bound. 

2. Upper bounds 

To get both upper bounds, we use the following trivial lemma. 
Lemma 2. Assume that (for some value of e) the powers ofT e satisfy 

Vn€N, HTfll < T(n), 

where the function T(n) is strictly decreasing, and r(n) n ^°°> 0. Then the dissipation time is bounded 
from above by 

where r* -1 ) is the inverse function of T . In particular, for the geometric decay T(n) = Cr n with 
r € (0, 1) ; C > 1, one obtains t* < nr^T + 



The upper bound in Eq. (|24j) comes from the obvious estimate 

< ||G e || n , 

on which we apply the lemma with C = 1, r = ||G e ||. 

To prove the second upper bound, we use the representation of T e n in terms of the resolvent: 



1 

2vri 



s r 



X n R t {X)d\ 



valid for any r > r sp (T e ). Thus for all r G (r sp (T e ), 1), one has 

< I |An|He(A)|||(iA| < sup \\R,(X)\\r n+1 = ^r" +1 . 
^TTjgr [A|=r d A r ) 

We then apply Lemma[21on the geometric decay for any radius r sp (T e ) < r < 1, with C = > 1. 

■ 

3.2. Consequences. We now use Theorem ^ in the case where U = Uf is the Koopman operator 
for some measure-preserving map F on T d with some specific dynamical properties. We recall jlUj 
that the map F is ergodic (resp. weakly-mixing) iff 1 is not an eigenvalue of Uf (resp. iff Up has 
no eigenvalue) on LqCF*). 

Proposition 2. For any measure-preserving map F and any noise generating function g, the 
dissipation time of T e diverges in the small-noise limit e — > 0. 

Proof. We skip the subscript F to alleviate notation. We only use the fact that U = Uf is an 
isometry. We shall prove by induction the following strong convergence of operators 

V/ G L 2 Q (T d ), Vn G N, ||T™/ - U n f\\^0. 

From Proposition this limit holds in the case n = 1. Let us assume it holds at the rank n — 1. 
Then we write 

T™/ = C/T™- 1 / + (G e - I)UT^- X f. 

From the inductive hypothesis, T™ _1 f ——*U n ~ l f , so that the first term on the RHS converges to 
U n f . Applying Proposition to the function U n f, we see that the second term vanishes in the 

limit e — > 0. From the isometry of U, we obtain that for any n > 0, ||T"|| e ~*°> l, so that r» £ ^°> oo. 
■ 

This "non-finiteness" of the noisy dissipation time allows us to prove another general result 
concerning the pseudospectrum of T t (this corollary is proved in the Appendix IA.3|) : 

Corollary 2. For any isometry U and noise generating function g, one has 

(27) 4(1)^0. 

This means that for any fixed 5 > 0, the pseudospectrum o~s(T e ) will intersect the unit circle for 
small enough e. 

In order to better control the growth of r* , we need more precise information on the noise and the 
dynamics. In the present section, we restrict ourselves to the dynamical property of weak-mixing. 

Corollary 3. Assume that the noise generating density g satisfies the estimates Q11JI or (\YS\i with 
exponent a G (0, 2] . If F is not weakly-mixing and at least one eigenfunction of Uf belongs to 
fji(Y d ) for some 7 > 0, then T e has slow dissipation time: 



-OA7) < 
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T*. 



Proof. Let h € H" l ("S d ) be a normalized eigenfunction of Up with eigenvalue A. Applying 
Proposition ^i), we get 

|| (A - T e )h\\ = \\{I - G e )h\\ < Ke> ha 

for some constant K > depending on g and h. This implies that ||i? e (A)|| > K€ \ac , therefore 
taking the supremum over |A| = 1 yields <i e (l) _1 > K ^,\ a ■ The lower bound in Theorem ^ then 
implies 

1 -e- 1 

(28) - — — < r*. 



Remark 1. Recall that if g satisfies with exponent a, then the dissipation time is also bounded 
from above, as shown in Corollary If one eigenfunction h has regularity H" 1 with 7 > a, then 
both Corollaries imply that the dissipation is simple, of exponent a. 

Remark 2. The above results can be stated in more general form: Uf does not need to be a 
Koopman operator associated with a map F . The result holds true for any isometric operator U on 
Lq with an eigenfunction of Sobolev regularity. 

The dependence of the lower bound in (|28[) on 7 can be intuitively explained as follows. In case of 
non- weakly-mixing maps the eigenfunctions of Uf are, in general, responsible for slowing down the 
dissipation. The rate of the dissipation is affected by the regularity of the smoothest eigenfunction. 
In principle, irregular functions undergo faster dissipation giving rise to slower asymptotics of r*. 
It is not clear, however, whether the actual asymptotics of the dissipation time will be slower than 
power law in case when all eigenfunctions of Uf on Lo(T d ) are outside any space H' y (T d ) with 
7 > 0. 

In Corollary 21 we have shown that for any map F and arbitrary small 5 > 0, the pseudospectrum 
cs(T e ) intersects the unit circle for sufficiently small e > 0. If F is not weakly-mixing, the spectral 
radius of T t (that is, the modulus of its largest eigenvalue) is believed to converge to 1 when e — > 0, 
and the associated eigenstate h e should converge to a "noiseless eigenstate" h. This "spectral 
stability" has been discussed for several cases in the continuous-time as well as for discrete-time 
maps on 

On the opposite, if F is an Anosov map on T 2 (see Section lo3|) . the spectrum of T e does not 
approach the unit circle, but stays away from it uniformly: r sp (T e ) is smaller than some ro < 1 for 
any e > [7j. Simultaneously, ||T e || — > 1, so we have here a clear manifestation of the nonnormality 
of T 6 for such a map. In some cases (see |25) and the linear examples of Section |SJ), the operator 
T e is even quasinilpotent, meaning that r sp (T e ) = for all e > 0. For such an Anosov map, the 
spectral radius of T e is therefore "unstable" or "discontinuous" in the limit e — > 0, while in the 
same limit the (radius of its) pseudospectrum o~s(T e ) (for 5 > fixed) is "stable". 

4. Universal lower bounds for the dissipation time of C 1 maps 

In this section we consider both noisy and coarse-grained evolutions. We start our discussion 
with general properties of the coarse-grained dissipation time. 

Proposition 3. Let F be a measure-preserving map. 

i) For arbitrary noise kernel, the coarse-grained dissipation time diverges as e — * 0. 

ii) If F is not weakly-mixing then r* = 00 for small enough e > 0. 

Proof. 

i) Similarly as in Proposition [21 we have 

\\T?f - U n f\\ = \\G e U3(Ge - I)f + (G e - I)U n f\\ < \\(G e - I)f\\ + ||(G e - I)U n f\\ -> 0. 
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ii) Let h G Lq(T^) be a normalized eigenfunction of Up, then 

\\ftoh\\ = \\G e U^(h + (G e -I)h)\\>\\G e h\\-\\G^(G € -I)h\\ 
> l-2\\(G e -I)h\\. 

Since the RHS above is independent of n, we see that ||rj n ^ || is close to 1 for all times and sufficiently 
small e > 0. ■ 

As opposed to the noisy case (see Prop.^), the coarse-grained evolution through a non-weakly- 
mixing map does not dissipate. Therefore there exists no general upper bound for f*. 

On the opposite, we will prove below a general lower bound for both coarse-grained and noisy 
evolutions, valid for any measure-preserving map F of regularity C . We note that Propositions |2] 
and |3t ) (which are valid independently of any regularity assumption) do not provide an explicit 
lower bound. 

First we introduce some notation. For any map F G C 1 , DF(x) is the tangent map of F at the 
point x G T d , mapping a tangent vector at a: to a tangent vector at F(x). Selecting the canonical 
(i.e. Cartesian) basis and metrics on T(T d ), this map can be represented as a d x d matrix. The 
metrics naturally yields a norm v G T x (T d ) i— > \v\ on the tangent space, and therefore a norm 
on this matrix: \DF(x)\ = max\ v \ = i \DF{x) ■ v\. We are now in position to define the maximal 
expansion rate of F: 

[i F = limsup ||DF n ||J/ n , where H-D-F^U = sup \{DF n ){x)\. 

n-*oo xe jd 

Since F preserves the Lebesgue measure, the Jacobian Jf(x) satisfies \Jf(x)\ > 1 at all points. In 
the Cartesian basis, Jf(x) = det(DF(x)), so that we have ||DF n (a;)|| > 1 for all x G T d , n > 0. 
One can actually prove the following: 

Remark 3. Although \(DF n )(x)\ and \\DF\loo may depend on the choice of the metrics, fip does 
not, and satisfies 1 < fj,p < ||I?-F||oo- 

From the definition of hf, for any [i> \xp there exists a constant A > 1 such that 

(29) Vn G N, WDF 11 ^ < Afi n . 

In some cases one may take /u = [Xf in the RHS. In case \xf — 1, ||-D-F n ||oo can sometimes grow as 
a power-law: 

(30) \\DF n \\ 00 < An^, n G N 

for some > 0, or even be uniformly bounded by a constant {(3 = 0). 

The relationship between, on one side, the local expansion of the map F and on the other side, 
the dissipation time, can be intuitively understood as follows. A lack of expansion (llDi 7 !^ = 1) 
results in the transformation of "soft" or "long- wavelength" oscillations into "soft oscillations", 
both being little affected by the noise operator G e . On the opposite, a locally strictly expansive 
map (H-DFlloo > 1) will quickly transform soft oscillations into "hard" or "short- wavelength" , the 
latter being much more damped by the noise. 

The following theorem precisely measures this relationship, in terms of lower bounds for the 
dissipation times. 

Theorem 2. Let F be a measure-preserving C 1 map on T d , and assume that the noise generating 
density g satisfies or M'J\) for some a G (0, 2]. 

i) If | IDF | |oo > 1, resp. fj,p > 1 3 then there exist a constant c, resp. constants [i > [if o,nd c, 
such that for small enough e, 

aAl , . i, _ a A 1 

(31) r* > - TTmFii \ ln(e ) + c, resp. r* > m(e ) + c. 

ln(||Z)F||oo) hi^x 
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If F is a C 1 diffeomorphism, then k31\) holds with \\DFWoo replaced by \\DFWoc A ^(F -1 ))!^, resp. 
with some fi > fip A fJ>p-i. 

ii) If HDi^lloo = 1 then T e has slow dissipation time, r* > g~( aA1 ). Jf the noise kernel satisfies 
the condition k!3]) for a G (0, 1], then the dissipation time is simple, r* ~ e _Q . 

Hi) If fip = 1 and \\DF n \\ OQ grows as a power-law as in Eq. ()3U|) with (3 > 0, then f* > e~( aA1 ^@ . 
if HD-F^lloo is uniformly bounded above by a constant, then f* = oo /or small enough e. 

Remark 4. JTiis theorem shows that classical systems on T d (i.e. C l diffeomorphisms) cannot have 
a dissipation time growing slower than Cln(e _1 ). In view of the results for toral automorphisms (cf. 
Proposition 4), this lower bound on the dissipation time is sharp and consistent with Kouchnirenko's 
upper bound on the entropy of the classical systems, namely all classical systems have a finite 
(possibly zero) Kolmogorov- Sinai entropy (Theorem 12.35. in PP, see also (5], |22| ). 

Proof of the Theorem. 

We will need the following trivial lemma (similar with Lemma EJ. 

Lemma 3. Assume that there exists some a > and a strictly increasing function 7(71), 7(0) = 
such that 

(32) Vn>l, ||T E n || > 1 - e Q 7 (ra). 
Then the dissipation time is bounded from below as: 

(33) T ,>J-1)(1ZJL V 

where 7 ( _1 ) is the inverse function 0/7. 

The same statement holds for the coarse-grained version. 

Our task is therefore to bound ||T e n || (resp. ||rj n ^||) from below. A simple computation shows 
that for any / G C°(T ), ||G £ /||oo < ll/lloo- Since convolution commutes with differentiation, for 
/ G C 1 we also have ||V((? e /)||oo < ||V/||oo- We use this fact to estimate the gradient of T € f: 

\MT e f)\\oo = \MG e U F f)\\oo 

< \\V(foF)\\ oc = \\(Vf)oF-DF\\ 00 

< ||(V/) o FWooWDFWn = IIV/IUIIDFIU. 
Repeating the above procedure m times, we get 

(34) ||V(T e m /)||oo < \\Vf\U\DF\\Z, ||V(C/ F 7r/)lloo < \\Vf\U\DF\\Z +1 . 

We now choose some arbitrary / G Cg(T d ), with ||/|| = 1. We first apply the triangle inequality: 

iiTj7ii = wGtUpTr'fw > wuprr'n - \\(G e - /^rrvii. 

To estimate the second term on the RHS we use the bound (|18|) and the estimate Q34|) to obtain 

> lirrVll - Ce^llv/iupFHSo. 

Applying the same procedure iteratively to the first term on the RHS, we finally get (remember 
11/11 = 1): 

n 

(35) \\T?\\ > ||TJ\f|| > 1 - C e QA1 ||V/|U \\DF\\Z- 

m=l 

The computations in the case of the coarse-grained operator are even simpler: 
\\fWf\\ = \\G e U%G e f\\ 

> 1 _ Ce aA1 ||V/|U - Ce aA1 ||V(a/)||oc||^ n ||oo 

(36) > l-2Ce aA1 ||V/|| 0C || J DF' 1 || 00 . 
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Notice that from the assumptions on /, ||V/||oo cannot be made arbitrary small, but is necessarily 
larger than some positive constant. We choose some arbitrary function, say / = with k = (1,0) 
which satisfies ||V/||oo = 2-7T. 

The estimate (jSEJ) has the form given in Lemma El The growth of the function r y(n) depends on 
whether H-Di^oo is equal to or larger than 1, which explains why the lower bounds are qualitatively 
different in the two cases. 

In case H-Di^Hco is strictly larger than 1, then the function j(n) grows like an exponential, 
therefore the lower bound is of the type 1)31)) . For the coarse-grained version, a growth of H-Di^oo 
of the type ()29]) yields the lower bound for f* in 1)3 1|). 

In the case H-Di^oo = 1, j(n) is a linear function, so that r* > gpr^i] — e~( aA1 ). 

In the coarse-grained version, if [ip = 1 and H-DF"!^ grows like in ()30|) with (3 > 0, the 
dissipation is slow: r* > Ce~( aA1 ^P . In the case where H-DF"!^ is uniformly bounded by some 
constant, the norm of the coarse-grained propagator stays larger than some positive constant for 
all times, so that for small enough noise f* is infinite. ■ 



5. An upper bound of the dissipation time for mixing maps 

For any two functions /, h € LqCT^), the dynamical correlation function for the map F is defined 
as the following function of n G N (see e.g. [3]): 

Cf,h(n) = Cj )fc (n) = m(fU$h) = {f,U$h) = {P$f,h). 

The same quantity may be defined for the noisy evolution: 

C e f , h (n) = m{fT?h). 

We recall that a map F is mixing iff for any /, h £ Lq, 

The correlation function can easily be measured in (numerical or real-life) experiments, so it is 
often used to characterize the dynamics of a system. 

To focus the attention, we will only be concerned with maps for which correlations decay in 
a precise way. We assume that there exist Holder exponents s*,s € M+, < s* < s together 
with some decreasing function T(n) = r StjS (n) with r(re) "~^°°> 0, such that for any observables 
/ £ Cg*(T d ), h € CQ(T d ) and for sufficiently small e > (sometimes only for e = 0), 

(37) VneN, \C} >h (n)\ < \\f\\c** \\h\\ c . T(n). 

In general, such a bound can be proven only if the map F has regularity C s+1 . The reason why we 
do not necessarily take the same norm for the functions / and h will be clear below. 
We will be mainly interested in the following two types of decay 

i) Power-law decay: there exists C > 0, [3 > such that, 

(38) r(n) = Cn-P. 

This behaviour is characteristic of intermittent maps, e.g. maps possessing one or several 
neutral orbits [1]. 

ii) Exponential decay: there exists C>0, 0<cx<l such that, 

(39) r(n) = Ca n . 

Such a behaviour was proved in the case of uniformly expanding or hyperbolic maps on the 
torus (see Section 0), as well as many more general cases jl]. 
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The central result of this section is a relationship between, on one side, the decay of noisy (resp. 
noiseless) correlations and on the other side, the small-noise behaviour of the noisy (resp. coarse- 
graining) dissipation time. The intuitive picture is similar to the one linking the local expansion 
rate to the dissipation: namely, a fast decay of correlations is generally due to the transition of 
"soft" into "hard" fluctuations of the observable through the evolution, which is itself induced 
by large expansion rates of the map. Still, as opposed to what we obtained in last Section, the 
following theorem and its corollary yields upper bounds for the dissipation time. 

Theorem 3. Let F be a volume preserving map on T d with correlations decaying as in Eq. \3T\) 
for some indices s, s* and decreasing function T(n), at least in the noiseless limit e = 0. Assume 
that the noise generating function g is ([s] + 1)- differentiate, and that all its derivatives of order 
\oc\x <[s\ + l satisfy 

\D a 9{x)\<T^M, 1*1 »!> 

with a power M > d. 

Then there exist constants C > 0, e Q > such that the coarse-grained propagator satisfies 



r(n) 

(/+*+*,. 



(40) Ve<e G , Vn > 0, ||T e (n) || < C 

e 

If the decay of correlations (j<-{7|) also holds for sufficiently small e > ( and assuming the Perron- 
Frobenius operator Pp is bounded in C s (T d )), then the noisy operator satisfies (for some constants 
C > 0, e Q > 0): 

(41) Ve< £o , Vn>0, ||T e l < C -*M-. 

From these estimates, we straightforwardly obtain the following bounds on both dissipation times 
(the assumptions on F and the noise generating function g are the same as in the Theorem): 

Corollary 4. I) If the correlation function satisfies the bound (|37|) for e = 0, then the coarse- 
grained dissipation time is well defined (f* < oo). Moreover, 

i) i/r(n) ~ n~P then there exists a constant C > such that 

Tif < Ce ~ 

ii) ifT{n) ~ a n then there exists a constant c such that 

~ d + s + s* 

t* < — n ; — m(e ) + c, 

I mcr l 

II) If Eq. (J37J) holds for sufficiently small e > 0, then 

i) ifT{n) ~ n _/3 ; there exists a constant C > such that 

d+s + s* 

r* < Ce — 

ii) ifT{n) ~ a n , there exists a constant c such that 

d+s + s* ! 

t* < — n i — m ( e ) + c - 

I m<T l 

Proof of the Theorem. 

1 st step: We represent the action of T" (resp. T^) on an observable / € Lq(T^) in terms of 
the correlation functions C e (n) (resp. C(n)). To do this we Fourier decompose both T™ +2 f and 
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fi = Upf, and use Eq. ©: 

T? +2 f = Yl {e h G e U F T^G e h)ej 

Q^j£Z d 

= E E fimGeej^FT^G^ej 

= E E fi(k)geU)Uk)(PFe j ,'lTe k )e j . 
o^jez d o^kez d 

(remember that g is a real function). A similar computation for the coarse-grained propagator 
yields: 

% n) f= E E /»^(j)ft(fe)(e i ,^e fe )e i . 
Taking the norms on both sides, we get in the noisy case: 

rr +2 /n 2 = E E hmPFej^e^umik) 2 

o=£jez 2 o^k& d 

^ E ( E l/i( fc )! 2 ) E \(PFej,T?e k )\ 2 \g e (j)Uk)\ 2 

(42) =^ ||T e ™+ 2 /|| 2 < H/ill 2 £ \Ch Fe _., ek (n)\ 2 \g(ej)g(ek)\\ 

0^j,k£Z d 

and in the coarse-graining case 

(43) ||TW/|| 2 <||/|| 2 Y, \C e ^, ek (n)\ 2 \g(ej)g(ek)\ 2 . 

These two expressions explicitly relate the dissipation with the correlation functions. 

2 nd step: We now apply the estimates (|37j) on correlations for the observables e k , e_j, Ppe^j. 
In the coarse-grained case, it yields (using simple bounds of the type of Eq. ) : 

Vj, k£Z d \{0}, |C e _.,e fc (n)| <C \j\ s \k\ s *T(n). 

In the noisy case, we need to assume that the Perron-Frobenius operator Pp is bounded in the 
space C s (T d ). This property is in general a prerequisite in the proof of estimates of the type (|57|l . 
so this assumption is quite natural here. 

Vj, k e Z d \ {0}, \Cp Fe _. >eh (n)\ < C \\P F e^\\o'\\e k \\ c ^T(n) 

(44) < C\\P F \\ C s\j\ s \k\ s *T(n). 

We insert these bounds on the decay of correlations in the expressions (|42H4H|) . for instance in the 
coarse-grained case we get: 

(45) Vn>0, ||T e (n) || 2 < CT(n) 2 (V( s+s *) ^ \ek\ s+s * g(ek)A . 

^ o^kez d ' 

3 rd step: We finally estimate the e-dependence of the RHS of the above inequality. Up to a factor 
e~ d , the sum in the brackets is a Riemann sum for the integral f \£\ s+s * g(£) 2 d£ < oo. A precise 
connection is given in the following lemma, proven in Appendix IA.4I 
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Lemma 4. Let f G C (M. d ) be symmetric w.r.t. the origin and decaying at infinity as |/(a;)| ;$ 
|a;| -A/ with M > d. Then the following estimate holds in the limit e — > 0: 

(46) e d V f{ekf = [ f + 0(e M ). 

kez* jRd 

Let m G N satisfy 2m < s + s* < 2m + 2 (notice that m < [s] since we assumed s* < s). From 
the obvious inequality 

Vx>0, x s+s * < x 2m + x 2m+2 , 

we may replace in the RHS of (|45[) the factor |efc| s+s * by |efc| 2m + |efc| 2m+2 . Applying Lemma 0] 
to the derivatives of g of order m and m + 1, we end up with the following upper bound, which 
proves the first part of the theorem: 

||fW|| 2 < Cr(n) 2 (^^^ d (^| 2m + |C| 2 ( m+1 ))5^) 2 ^ + 0( e A - ^ 
< C> F(n)2 Hall 4 R 

The computations follow identically for the case of the noisy operator, yielding the second part of 
the theorem. ■ 

6. Some examples of mixing maps 

In this section we apply the results of the last two sections to several classes of Lebesgue measure 
preserving maps which have been proven to be mixing, with various types of correlation decays. 

Remark: In general, an expanding or hyperbolic map on the torus does not preserve the 
Lebesgue measure, so the first step is to precise with respect to which invariant measure one wants 
to study the ergodic properties. In the "nice" cases, one can prove the existence and uniqueness 
of a "physical measure", which is ergodic for the map F, and then study the (noisy or noiseless) 
mixing properties with respect to this measure. As pointed out in the introduction, in the present 
paper we only consider maps for which the physical measure is the Lebesgue measure. 

6.1. Decays of correlations. We briefly summarize the results concerning mixing maps, mostly 
relying on the review ,4j and the book [Sj. 

Let us first consider the noiseless correlations. A common route to prove that a map F is 
exponentially mixing consists in identifying an invariant space B of densities, on which the Perron- 
Frobenius has the following spectral structure: except for the eigenvalue 1 associated with the 
constant density, the rest of the spectrum (on the subspace Bo of zero-mean densities) is contained 
inside a disk of radius er < 1 centered at the origin. More precisely, the spectrum on Bo may consist 
in isolated eigenvalues (called resonances) {A,} and of some essential spectrum in the disk of radius 
r ess < a. Pp is said to be quasicompact on B. This spectral structure implies that there is some 
constant C > such that for any / G Bo, h G Bq, 

(47) Vn > 0, |C /lfc (n)| = \{h,P$f) B . )B \ < \\h\\ B * ||P£|| Bo ||/|| B < C \\h\\ B * \\f\\ B ° n ■ 

For the maps we study, the spectrum of Pp on Lq intersects the unit circle, so B cannot simply 
be the Hilbert space L 2 . Depending on the properties of the map, B can be a Frechet space of 
analytic functions, a Banach space of bounded variation, Holder or C s functions (see Section f(i.2|) : 
it may also be a space of generalized functions lying outside L 2 (see Section ESI • No matter how 
complicated B can be, in general there exist Holder exponents < s* < s such that C s (resp. 
C s *) embeds continuously in B (resp. in its dual B*). As a result, the upper bound (|47|) can be 
specialized to functions / G Cq, h G Cq* as follows: 

(48) Vn>0, IC^^C^)! < C H^llo-* ||/||C7- 
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This is the form of upper bound we used in Theorem |21 

This strategy of proof has been applied to several types of maps, including the (noninvertible) 
expanding maps and the Anosov or Axiom-A diffeomorphisms on a compact manifold. Exponential 
decay of correlations has also been proven (using various methods) for piecewise expanding maps on 
the interval, some nonuniformly hyperbolic/expanding maps, some expanding or hyperbolic maps 
with singularities. 

Other types of decay occur as well: for instance, a polynomial decay of correlations Cfh(n) < 
was shown to be optimal for some "intermittent" systems, like a one-dimensional map expanding 
everywhere except at a fixed "neutral" point (such maps are sometimes called "almost expanding" 
or "almost hyperbolic"). 

There exist fewer results on the decay of correlations for stochastic perturbations of deterministic 
maps, like our noisy evolution T e . In general, one wants to prove strong stochastic stability, that 
is stability of the invariant measure and of the rate of decay of the correlations in the small-noise 
limit. In our case, only the second point needs to be proven, since the Lebesgue measure remains 
invariant after switching on the noise. 

Stochastic stability has been proven for smooth uniformly expanding maps |Hj (see next sub- 
section), some nonuniformly expanding or piecewise expanding maps. It has been shown also for 
uniformly hyperbolic (Anosov) maps on the 2-dimensional torus (see subection l6.3|) . In all those 
cases, the mixing is exponential, so the stability of the decay (|48|) means that for small enough 

e > 0, there exists a radius a t e ^°> a such that for any / G Cq, h € Cq* , 

(49) Vn > 0, |<%(n)| < C \\h\\ C s* \\f\\c° ^. 

Here the constant C > can be taken independent of e. 

In the next two sections, we describe in more detail the cases of smooth uniformly expanding 
maps and Anosov diffeomorphisms on the torus. 

6.2. Smooth uniformly expanding maps. Let F be a C s+1 map on T d (with s > 0). Assume 
there exists A > 1 such that for any x £ T d and any v in the tangent space T x T d , \\DF(x)v\\ > X\\v\\ 
(we assume that A is the largest such constant). Such a map is called uniformly expanding. In 
general, it admits a unique absolutely continuous invariant probability measure; here we restrict 
ourselves to maps for which this measure is the Lebesgue measure. 

Ruelle (201 proved that the Perron-Frobenius operator Pp of such a map is quasicompact on 
the space C s (T d ), and that its essential spectrum is contained inside the disk of radius r\ = A~ s . 
In general, one has little information on the possible discrete spectrum outside this disk (upper 
bounds on the decay rate have been obtained in the case of an expanding map of regularity C 1+v 
[2]). Strong stochastic stability for such maps was proven in [2], with a more general definition of 
the noise than the one we gave. 

For all these cases, one can take = 0, since the continuous functions are continuously embedded 
in any space (C s )*. 

Case of a linear expanding map. We describe the simplest example possible for such a map, 
namely the angle-doubling map on T 1 defined as F(x) = 2x mod 1. This map is real analytic, with 
uniform expansion rate A = 2. Due to its linearity, the dynamics of this map (as well as its noisy 
version) is simple to express in the basis of Fourier modes e k (x) = e 2mkx -_ 

Vfc € Z, U F e k = e 2k 

=> T e e k = g{ek)e 2k 

n 

^T?e k = [Hg(e^k)]e 2 n k 
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The computation is even simpler for the coarse-grained propagator: 

fV*e h = g(k)g(2 n k)e 2 n k . 

To fix the ideas, we consider the a-stable noise <?(£) = e - '^'" for some < a < 2. One easily checks 
that for any n > 1, 



e a - 



||T e n || = ||T e n ei|| =exp 

IIT^H = exp{-e°(2 na + 1)}. 

For any e > 0, these decays are super-exponential: the spectrum of T e on L\ is reduced to {0} for 
any e > (the spectrum of Up is the full unit disk). From this explicit expression, we get both 
dissipation times: 

(50) n = j-L Me" 1 ) + 0(1), f* = ^ ln^ 1 ) + 0(1). 

For this linear map, H-Di^Hoo = = 2, so this estimate is in agreement with the lower bounds (|31j) . 
the latter being sharp if a £ [1, 2]. On the other hand, In 2 is also equal with the Kolmogorov-Sinai 
(K-S) entropy h{F) of F. Therefore, for this linear map the dissipation rate constant exists, and 
is equal to j^p^j- 

To compare these exact asymptotics with the upper bounds of Corollary we estimate the 
correlation functions Cj^in) on the spaces C^T 1 ). We give below a short proof in the case s > \. 
We will use the following Fourier estimates |33| : 



aoo, ,V/ e c^t 1 ), vmo, |/(fc)|<c 

Therefore, writing the correlation function as a Fourier series, we get: 



\ Pm ? = E \h^)\ 2 < E ( c w§ 



(si) =MWII < C { , 

This estimate yields a decay of the correlation function as in Eq. I|48|l. with a rate a = 2~ s and 
s* = 0. One can check that this rate is sharp for functions in C s : indeed, any z € C, \z\ < 2~ s 
is an eigenvalue of Pp on that space. Applying the Corollary0J/ ii), we get that for any s > 1/2, 
there exists a constant c such that 

(52) f* < i-ti ln^j + c 

sm2 

for sufficiently small e. The exact dissipation rate constant l/ln2 is recovered only for large s. 

A straightforward computation shows that the estimate ()51|) also holds if one replaces Pp by 
Pp o G e ; hence the noisy correlation function dynamics satisfies the same uniform upper bound 
as the noiseless one, with the decay rate a e = 2~ s . As a result, the upper bound on r* given by 
Corollary ii) is the same as in Eq. (|52JI . 



6.3. Anosov diffeomorphisms on the torus. We recall that a diffeomorphism F : T d i— ► T d is 

called Anosov if it is uniformly hyperbolic: there exist constants A > and < A s < 1 < X u such 
that at each x £ T d the tangent space T x T d admits the direct sum decomposition T x T d = E s x ®El 
into stable and unstable subspaces such that for every n € N, 



(D X F)(E X ) = E Fx , \\(D x F n ) lEg \\ < AX 
{D X F){ED = E Fx , \\(D x F- n ) lE J\ < AX- 
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These inequalities have obvious consequences on the expansion rates of F and F , for instance they 
imply llDi 7 !!^, > H-DF"!^ > A _1 A™. As a consequence, the quantities of interest in Theorem [3 i) 
satisfy 

IIDFIIooAHDF- 1 !^ > a u aa;\ 

/ijr A fXp-l > A^AA^ 1 - 

All these expansion rates are > 1, so both noisy and coarse-grained dissipation times admit loga- 
rithmic lower bounds as in Eq. (|31[) . 

Exponential mixing has been proven for Anosov diffeomorphisms of regularity C 1+v (0 < rj < 1) 
by Bowen [S], using symbolic dynamics; the exponential decay is then valid for Holder observables 
m C for some < ?/ < rj. Because we are also interested in the noisy dynamics, we will refer 
to a more recent work j2] concerning C 3 Anosov maps on T d , which bypasses symbolic dynamics. 
The authors construct an invariant Banach space B of generalized functions on the phase space, 
such that the Perron-Frobenius operator is quasicompact on this space. One subtlety (compared 
with the case of expanding maps) is that B explicitly depends on the (un)stable foliations of the 
map F on T d . Vaguely speaking, the elements of B are "smooth" along this unstable direction 
but can be singular ("dual of smooth") along the stable foliation E^,. The space B is the 
completion of C 1 (T c( ) with respect to a norm || • \\q adapted to these foliations. This norm, and 
therefore B, are defined in terms of a parameter < [3 < 1, the choice of which depends on the 
regularity of the unstable foliation. In general, the latter is r-Holder continuous, for some exponent 
< r < 2. Then, the authors prove that if one takes (3 < r A 1, then the essential spectrum of 
Pf on B has a radius smaller than rp = max(A~ 1 ,Af). This upper bound is sharper if (3 can be 
taken close to 1, that is, if the foliation is C . This is the case for smooth area-preserving Anosov 
maps on T 2 , for which the foliations have regularity C 2 ~ s for any 5 > The operator Pp may 
have isolated eigenvalues (resonances) 1 > |Aj| > rp, corresponding to eigenstates in Bq which are 
genuine distributions Lq. There is (to our knowledge) no simple general upper bound for the 
largest resonance |Ai| in terms of the expansion parameters (A u , X s ). 

By construction, the space C 1 (T d ) embeds continuously in both B and its dual B*, so that one 
can take s = = 1 in Eq. (|48|) . Therefore, for any erg > max(|Ai|, rp), there is some constant 
C > such that for any f,he C£(T d ), 

(53) Vn>0, \C f , h (n)\<C\\h\\ cl \\f\\ c ,a n p. 

In the proof of Theorem 01 (Step 3), for the case s = s* = 1 we only need to assume that the noise 
generating function g is C 1 with fast-decaying first derivatives. The fast decay implies that the 
first moment of g is finite (that is, one can take a > 1). 

The noisy propagator G € Pf is also analyzed in j2j ■ If the unstable foliation has regularity C 1+T> 
with 77 > (for instance for any C 3 Anosov diffeomorphism on T 2 ), and assuming that the noise 
generating function g € C 2 (M. d ) has compact support 1 , the authors prove the strong spectral stability 
of the Perron-Frobenius operator Pp on any space B defined with a parameter (3' < rj. Notice that 
this constraint on /3 is stronger than in the noiseless case, where one could take any f3 < 1; the 
spectral radius a pi may accordingly be larger than a p. Modulo the replacement of ap by ap/, the 
estimate (|53|) therefore applies to the noisy correlation function Cj h (n) as long as e is small enough. 

We collect below the results regarding the dissipation time of C 3 Anosov maps on the torus. 

Theorem 4. Let F be a volume preserving C 3 Anosov diffeomorphisms on T , and let the noise 
generating function be C 1 with fast decay at infinity. 



The condition of compact support could probably be relaxed to one of fast decrease at infinity (C. Liverani, 
private communication) 
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I) Then there exist fi>\ u /\\ s 1 ,0<a<l and C > such that the dissipation time of the 
coarse-grained dynamics satisfies 

^- He- 1 ) - C < f* < ln ( e ~ x ) + c, 

liifi | In cr | 

II) If in addition F has C l+V -regular foliations, and g £ C 2 (M. d ) is compactly supported, then 
there exist a < a < 1 and C such that the dissipation time of the noisy dynamics satisfies 

He" 1 ) - C < r* < ^4 m ( e_1 ) + C 



lnH-D-Flloo |l ncr l 

6.4. Ergodic linear automorphisms of the torus. In this section we describe examples of 
Anosov maps for which the dissipation time can be precisely determined. One can even compute 
the dissipation rate constant; we will mention the connection of the latter with the Kolmogorov- 
Sinai entropy. After the simple example of a (generalized) cat map on the 2-dimensional torus pQ, 
we recall the results obtained in for d-dimensional ergodic automorphisms. 

Throughout this section, we take a noise generating function of the type (fTl)) for a certain 
a € (0,2], with Q = I: 

(54) g eia (x) := £ e'W* e k {x). 

kez d 

Notice that this function is not compactly supported, therefore it does not satisfy stricto sensu 
the assumptions required in [7j to prove the strong spectral stability of Pp (see the footnote after 
Theorem 0} . 

2-dimensional cat map. We take A € SX(2,Z) with |TrA[ > 2 and define the dynamics as F(x) = 
Fa(x) = A l x mod 1 (A* is the transposed matrix of A). This map is of Anosov type, and preserves 
the Lebesgue measure. The dynamics is easy to express on the Fourier modes: 

Up e k = eAk 

T?e k = e-HUM^^ 

Similarly, 

\\fW\\ = exp (-e a min (\k\ a + \A n k\ a )] . 

Let us call A, A -1 the eigenvalues of A, with the convention |A| > 1. One can easily show that there 
are two constants < C\ < C2 such that for any n > 0, 

n 

Ci \X\ na/2 < min V \A l k\ a < C 2 \\\ na/2 , 

O^fcGZ 2 ^ 

d \\\ na/2 < min (|fc|° + \A n k\ a ) < C 2 \X\ na/2 . 
The above estimate yields the following asymptotics in the small-e limit: 

2 Me" 1 ) + 0(1), f„ = -4-ln(e~ 1 ) + 0(l). 



* ln|A| v ' w ' * ln|A| 

As in the case of linear expanding maps, the dissipation rate constant exists, and seems related 
with the K-S entropy of the linear map h(F) = In |A|. 

Let us compare these exact asymptotics with the bounds obtained in the previous sections. 
The cat map being linear, one has H-D-FHoo = \\A t \\ > |A|, H-DF" 1 ^ = > |A| (with 

equality iff A is symmetric). Since A is diagonalizable, we have ||DF n || 00 = ||(A') n || ~ |A|", so that 
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Hf = Hf- 1 = \M- Therefore, the lower bounds for the dissipation times given in Theorem^J are 
strictly smaller than the exact rates derived above, they differ from the latter by a factor < 1/2. 

To estimate the dynamical correlations, we use Fourier decomposition and proceed as in 
we construct a "primitive subset" S of the Fourier plane Z 2 \ {0}, such that each Fourier orbit 
intersects S once and only once. This subset looks like the union of 4 angular sectors of the plane, 
and has the following crucial property: 

(55) 3c>0, VfcGcS, VneZ, \A n k\ > c|A| |n| \k\. 

The correlation function between two observables f,hE Lo(T d ) may be written as: 

Cf,h(n)= f(-k)HA n k) = ^2^2f(-A l k)h(A l+n k). 

(Vfcez 2 kes lez 

Let us assume s > 1, and take / £ Cq(T 2 ): its Fourier coefficients decrease for all ^ k G 7? as 



|/(fc)l<c^" |r " 



\k\ s 

Using this decrease as well as the estimate (|55|). we find that 



\C f , h (n)\<C 2 \\f\\ C s \\h\\c-2^X. 



The sum over k converges, and the sum over I is bounded from above by ttou- This implies that for 

any s > 1, the correlation function decreases as in Eq. Q48|) . with the rate a = |A| _S . Actually, this 
decrease holds as well in the case < s < 1, but the proof is different. From there, Corollary 0]f ii) 
yields the upper bound 

(56) f. < He' 1 ) + c. 

This method can be straightforwardly adapted to prove that the noisy correlation function 
Cf h (n) decreases as fast as the noiseless one (in case s > 1). As a result, we obtain the same upper 
bound for as for f*. We notice that, as in the case of the angle-doubling map, the constant in the 
upper bound converges to the exact rate only in the limit s ^> 1, while for s = 1 (cf. the discussion 
preceding Thm.Q, the constant is twice larger. 

Toral automorphisms in higher dimensions. We consider toral automorphisms F = F& given by 
matrices A € SL(d,Z); such an automorphism is ergodic iff none of the eigenvalues of A are roots 
of unity. The K-S entropy h(F) of F is given by the formula jSJ 

(57) h(F) = m l A il> 

|Aj|>l 

where Xj are the eigenvalues of F. All ergodic toral automorphisms have positive entropy. We 
denote by P the characteristic polynomial of the matrix A and by {Pi, ...,P r } the complete set of 
its distinct irreducible factors (over Q). Let dj denote the degree of the polynomial Pj and hj the 
K-S entropy of a toral automorphism with the characteristic polynomial Pj. For each Pj we define 
its dimensionally averaged K-S entropy as 

(58) hj = h f. 

Uj 

For the whole matrix F we define its minimal dimensionally averaged entropy (denoted h(F)) as 

h(F) = min hj 
3=1,— >r 
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Notice that for an irreducible matrix A, that is a matrix admitting no proper invariant rational 
subspace, this quantity reduces to ^p-- 

Our notation regarding the noise parameter slightly differs from the one used in |14j (one has to 
replace e 2a by e); the results obtained there read as follows: 

Proposition 4. ^3] Let F = Fa be a toral automorphism on T d , and assume the noise kernel to 
be as in Eq. \54\ l- Then, in the limit of small noise, 

i) both dissipation times have a power-law behaviour iff F is not ergodic. 

ii) both dissipation times have logarithmic behaviour iff F is ergodic. 
Hi ) if F is ergodic and A is diagonalizable then 

r* ss f* m - ln(e _1 ). 

h(F) 

We end this section with a remark that the small-noise asymptotics of the dissipation time is 
insensitive to a super-exponential decay of the correlation functions. We illustrate this fact by the 
following result on the decay of correlations for d-dimensional toral automorphisms, which can be 
proven along the same lines as the above proposition: 

Proposition 5. Let F be a diagonalizable ergodic toral automorphism and A any constant such 
that < A < h(F). Then for any f,h£ Lg(T M ) 

CUn)< H/ll \\h\\e~^ n 

Let f,h£ G e (LQ(T 2d )) be smooth observables. Then 

Cf^KWG-'fUG^hWe-^. 

7. Conclusion 

We have investigated the effect of noise or coarse-graining on the dynamics generated by a con- 
servative map, in particular the connection between the speed of dissipation of the noisy dynamics 
and the spectral and dynamical properties of the underlying map. 

We restricted ourselves to the case of volume preserving maps on the d-dimensional torus. The 
choice of the torus allowed us to use the Fourier transformation in our proofs, that is, harmonic 
analysis on that manifold. Most of our results can certainly be generalized to volume preserv- 
ing maps on other compact Riemannian manifolds. The choice of the torus was also guided by 
the existence of simple volume preserving Anosov maps on it, most notably the linear examples 
presented in Sections 16.21 16.41 As explained at the beginning of Section |H1 one may also want to 
extend the results to maps which do not leave invariant the Lebesgue measure, but still admit a 
'physical measure', as is the case for uniformly expanding or Anosov maps (the physical measure is 
of SRB type, that is, its projection along the unstable manifold is absolutely continuous jS]). The 
noisy perturbations of these maps have been analyzed as well Ej , in particular the strong spec- 
tral stability of the Perron-Frobenius operator holds under the same conditions as for the volume 
preserving maps. Although the equilibrium measure is more complicated, it might be possible to 
define and study a dissipation time in this more general framework. 

As we explained in Corollary 13 the dissipation of a non-weakly-mixing map is governed only 
by the nontrivial eigenstates of the Koopman operator, the speed of dissipation depending on the 
smoothness of these eigenstates. The asymptotics of the noisy dissipation time is generally a power 
law in e , while there remains a possibility of faster dissipation if all eigenstates are singular 
enough (see Remarks after Corollary ^ . 

The more interesting results concern the mixing dynamics, in particular when the mixing occurs 
exponentially fast. We then proved that both noisy and coarse-graining dissipation times behaved 
as the logarithm of e -1 in the small noise limit. This dependence can be understood through the 
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time evolution in Fourier space: a mixing map transforms long wavelength fluctuations into short 
wavelength ones, the latter being damped fast by the noise operator. This evolution in Fourier space 
is typical of uniformly expanding/hyperbolic maps; it is already responsible for the fast decay of 
dynamical correlations. The link between this decay and the dissipation is explicited in Theorem|31 
and its application to Anosov maps is given in Theorem |IJ 

In this context, we were unable to solve the problem of the existence and value of the dissipation 
rate constant (that is, the prefactor in front of ln(e -1 )). We obtained lower, resp. upper bounds 
for this constant, in terms of the local expanding rate, resp. the rate of decay of correlations. It is 
not clear whether this constant can in general be related with the measure-theoretic (KS) entropy 
of the map, as is the case for linear automorphisms (modulo some algebraic subtleties 1141). 

Although we used the spectral estimates of Theorem ^mostly for the case of non- weakly-mixing 
maps, it also makes sense to use that theorem in the reverse direction in the case of mixing maps, 
that is, deduce (pseudo)spectral properties of the noisy propagators, starting from the dissipation 
time estimates obtained in Corollary The analysis of the pseudospectrum of T t for mixing maps 
would complement the spectral one [ZIES]- We did not enter this aspect in the main text, because 
our attention was devoted to obtaining information on the dissipation time. 

The results of the present paper concern classical (i.e. non-quantum) dynamical systems. The 
quantization of both linear and nonlinear maps on a symplectic (even-dimensional) torus as a phase 
space has been studied in a number of works ^Hl El EH [201 EES E2] • Several recent studies deal with 
some form of noise, or decoherence, in discrete-time quantum dynamics |9l I15 | 124" ! I16 | 125] . While 
most of these works concentrate on spectral or entropic properties of noiseless/noisy dynamics, the 
long time behaviour of the quantum system can also be studied from the dissipation time point of 
view. The quantum setting provides a natural framework for extension of the present work, which 
we will address in a separate paper [13] . 

Acknowledgments: This project started during a stay of S.N. at the Mathematical Sciences 
Research Institute (Berkeley), the support of which is gratefully acknowledged. The authors would 
like to thank Prof. S. De Bievre and Prof. B. Nachtergaele for putting them together in contact. 



A.l. Proof of LemmaEJ We use the following upper bound: for any a £ (0, 2], there is a constant 
C a such that 



Appendix A. Proofs of some elementary facts 



(59) 



ViGR 



< 1 -cos(2vrx) < C a \x\ a . 

2 

cos(x) ss Tj- for small x. We simply apply these estimates to 



Besides, one has the asymptotics 1 
the following integral: 




In the case g admits a second moment, we have in the limit £ 



0: 




where we have used the notation £ 



l|y for any £ / 0. 
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A. 2. Proof of Proposition ^ . The statement i) is standard in the context of distributions |3L)l 
p. 157]. In our case, assume that / G L 2 is normalized to unity and consider an arbitrary small 
5 > 0. Since / E L 2 (T d ), there exists K > s.t. ^|fc|>_ft" l/COl 2 < ^- Since g is continuous and 
<?(0) = 1, there exists rj such that (1 — g(£)) 2 < 6 if |£| < 77. Thus using spectral decomposition © 
of G e , we obtain for all e < -ft 

(60) i|G e / - /f = e a - 5(^)) 2 i/(fc)i 2 < 5 E i/»i 2 + E i/»f 2 ^ M - 

fcGZ d fc|<if fc|>ii" 

To prove the next statement, first notice that if g satisfies the estimate (|12[) for the exponent 
a, it also satisfies it for the exponent 7 A a. Using once again spectral decomposition of G e , and 
applying the estimate (|12[) with the latter exponent we get 

na/-/n 2 < E( c 7AaM 7Aa |6fer a ) 2 i/>)i 2 

fcez d 

(61) < (C, Aa M jAa ) 2 e 2 ^ E |fc| 2(7Aa) l/>)! 2 

fcez d 

< (C 7Aa M 7Aa ) 2 e 2 ^||/||l, 7A(l . 

To obtain the last statement, we notice that any / € C 1 (T c( ) is automatically in H l (T d ), and 
that its gradient satisfies 

\\vf\\i > iiv/11 2 = 4vr 2 e ifciv»i 2 > ^ 2 e ifci 2(1Ao) i/»i 2 . 

k£Z d kez d 
The inequality (|oT|) with 7 = 1 then yields the desired result. ■ 

A. 3. Proof of Corollary l2l We prove the limit d e (l) £ ~*°> by contradiction. Assume that there 
is some constant a G (0, 1) such that for all e > 0, the distance d e (l) > a. We will show that the 
following triangle inequality holds: 

(62) Ve > 0, (4(1 - a/2) > a/2. 

First of all, notice that the assumption d e (l) > a means that for any A £ S , ||i? e (A)|| < a . We 
apply the following identity [3*0] : 

R e (\>) = R e (X){l + E(A - A') n ^(A)"} 

n>l 

with A' = rA, for 1 — a < r < 1. Taking norm of both sides yields the bound ||i? e (A')|| < r _^_ a ^ , 
uniformly w.r.t e. Since this upper bound holds for any |A'| = r, it shows that the spectral radius 
r sp (T e ) < 1 — a, and proves (|6*2*|) by taking r = 1 — a/2. We can now use ()62p in the upper bound 
(|25|) of Theorem^ this e-independent upper bound shows that remains finite in the limit e — ► 0, 
which contradicts Proposition [2j ■ 

A. 4. Proof of Lemma |1J Considering its decay at infinity, the function / is automatically in 
L 2 (M. d ). The function f 2 is the Fourier transform of the self-convolution f * f. Therefore, using 
the parity of / and applying the Poisson summation formula to the LHS of (|46|) yields 



(63) e d E f(ek) 2 = [ / 2 (£K + E (/*/)(" 



f*f)\ 

-M 



A simple computation shows that (f*f)(x) also decays as fast as |x| — . This piece of information 
is now sufficient to control the RHS of (|63j) . yielding the result, Eq. (|46|) , ■ 
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A. 5. Proof of Corollary [SJ It was showed in |14j that for any 5 > in case of Gaussian noise 

(64) \\T?\\<e-* a ' 3(1 - g)HF)n . 

Using this estimate one immediately gets 

q» = (f,T?g) < ||/|||b||||rf|| < ||/|||b||e- e2A ". 

Now let / = G t fo and g = G e go. Since the estimate (|6"4"]l holds also in coarse-grained version, we 
have 

C f>g (n) = (f,U£g) = (Gjo,UZG e g ) = (f ,f^g ) 



2 \ n 



< !l/ollll5olll|T^||<||/olllbol|e^ A . 
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